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Abstract 



Partial wave amplitudes have been extracted from accurate data on pp — ► 
7r~7r + , combined with earlier data on pp — > 7r°7r°, by a method which in- 
corporates the theoretical constraints of analyticity and crossing symmetry. 



^ ■ The resulting solution gives a good fit to the annihilation data and is also 

consistent with the wealth of information in the crossed channel irN — ► 7riV. 

The partial wave amplitudes show evidence for resonances in all partial waves 

with J < 5, at least one of which, a J = + state, (and possibly another with 

J = 1~) is unlikely to have a simple qq structure. 
13.75.Cs, 11.80.Et, 14.40.Cs, 25.43.+t 
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I. INTRODUCTION 



There have been a number of attempts to obtain information on meson resonances from 
data on the reaction NN — > -kit. Earlier work was restricted to testing the compatibility 
of data with assumed combinations of Breit-Wigner resonances and smooth backgrounds. 
These analyses used data only on the channel 7r~7r + and the results, in the main, are in 
strong disagreement with later tt°tc data [[| . An alternative approach is to use the Barrelet 
zero method Q to translate the dips observed in the angular distributions (cross-sections 
and asymmetries) into complex zeros of the scattering amplitudes. This method requires a 
number of assumptions to obtain smooth zero trajectories and to fix the overall phase of the 
solution. The latter is crucial, because in the absence of an optical theorem the method only 
determines the moduli and relative phases of the amplitudes at each energy. One choice is 
to assume a Breit-Wigner form for a particular partial wave, but different choices of input 
phases, however plausible, lead to different solutions for the partial waves, as can be seen 
by comparing the results of references || and ||. 

An alternative method of obtaining amplitudes which overcomes the objections to earlier 
work, is to exploit analyticity via dispersion relations and determine the overall phase by us- 
ing crossing symmetry to relate the annihilation data to known ttN — > ttN elastic scattering 
amplitudes. This in principle overcomes the problem of the absence of an optical theorem 
and ensures that the resulting amplitudes not only fit the annihilation data, but are also 
consistent with the wealth of information on the crossed channel, ttN — > ttN. Dispersion 
relations at fixed-t or fixed-w are unsuitable because of the need to have information in un- 
physical regions, but this can be avoided by writing dispersion relations along hyperbolae in 
the Mandelstam plane [[jj . In a previous paper [§] we presented a largely model-independent 
set of invariant amplitudes, which by construction satisfy analyticity and crossing symmetry 
and simultaneously give an excellent fit to all pp — > hit data (both in charged and neutral 
channels) existing at the time HHHU and to amplitudes for ttN — > ttN. Subsequently [H| 
we showed that the structure in the invariant amplitudes, which reflects that seen directly in 
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the experimental data, was due to the existence of resonances in the partial wave amplitudes 
and we presented evidence for resonance activity in four states with spins 1, 2, 3 and 4. 

Since our earlier work, new differential cross-section and asymmetry data have been 
obtained on the 7r~7r + channel at 20 center-of-mass energies in the range 1.91 MeV to 2.27 



MeV [IT] using the LEAR facility at CERN. These data are consistent with the earlier data 
H, but are more accurate and also extend asymmetry measurements to lower momenta 
than previous experiments. However, the new experiments cover a slightly smaller range of 
energy. Here we present an analysis of the new data, plus the older 7r°7r° data 101, usine; 
as constraints the invariant amplitudes obtained in ||, thus indirectly imposing analyticity 
and crossing symmetry on the solution. We compare our results with other recent analyses 



12-14 at the end of this paper. 



II. PARTIAL- WAVE AMPLITUDES 

The analysis of || was performed in terms of the usual ttN invariant amplitudes A and 
B. In the t-channel, NN — > tttt, it is more convenient to work with B and 

C = -A + m(^\ cos 9 t B, (1) 

where q(p) is the center-of-mass momentum in the im(NN) channel; M is the nucleon mass; 
and 9 t is the scattering angle in the t-channel. In the helicity basis, these two invariant 
amplitudes may be expressed in terms of partial-wave amplitudes by the expansions 

B\t, cos 6 t ) = —Y 2J+1 P' (cos 9 t ) Fj (t) (2) 

and 

C 1 (t, cos 9 t ) = ^ Y, ( 2 J + !) P -J ( cos °t) N J (*) • ( 3 ) 
P J 

In these relations, Nj and Fj are the partial-wave helicity amplitudes for definite t-channel 
isospin / corresponding to t-channel helicity non-flip and flip, respectively, and t is the square 
of the center-of-mass energy in the t-channel. By Bose statistics, 



N}(t) = F}(t) =0 for J even (4) 

and 

Nj(t) = F°j(t) = for J odd. (5) 

In addition, F® = 0. The helicity amplitudes are normalized such that the integrated cross- 
section for NN — > 7T7T for a definite isospin I is 



- J = 2- (^)E(2^+1){K + |FJ 



(6) 



An alternative, but equivalent, way of expressing the partial- wave content of the invariant 
amplitudes is to use the LS basis. If we label the new amplitudes by JL, i.e. Hjl, then, for 
a fixed value of isospin J, they are related to Nj and Fj by 

Hj,l=j-i = (^t) V2 {^J + 1Fj ~ ^JNj) (7) 



and 



Hj, L=J+ i = j) {^/JFj + VJTlNj} (8) 



with 



^ = 7r (^)E(2^ + i){l^./+l 2 + |i/ J -i 2 }, 



(9) 



where Hj± = H JL=J ± 1 . We will use both sets of amplitudes in the following discussion. 

In ||, the invariant amplitudes were obtained along families of hyperbolic curves in the 
Mandelstam plane defined by a parameter 

(M*-m 2 ) 2 - 3 (S- 3 -t) 
^~ t-4m 2 { } 

where m is the pion mass; s is the square of the center-of-mass energy in the s-channel 

ttN — > ttN; and S = 2 (M 2 + m 2 ). Values of £ were chosen so that the hyperbolae covered 

the t-channel region being analysed while staying almost entirely within the physical s- 

channel. To obtain amplitudes at the experimental energies, those from ref. M were linearly 

interpolated in t along the hyperbolae at fixed-£. 



III. FITS TO DATA AND AMPLITUDES 



The data fitted consisted of experimental differential cross-sections and asymmetries for 



pp — > 7r~7r + at 20 momenta in the range 360 MeV/c to 1550 MeV/c [I I] and, at each 
momentum, the invariant amplitudes from ref. |J at a grid of ^-values (or equivalently, a 
grid of values of cos#t). In addition, differential cross-sections for pp — ► 7r°7r were included 
by interpolating the measured data of ref. to the momenta of the charged channel data. 
At each momentum we minimized the sum of the values for x 2 per data point for each type of 
data. The parameters are, at each t- value, the complex amplitudes Nj and Fj. The number 
of amplitudes used in the partial-wave expansions is dictated by the form of the invariant 
amplitudes and the need to ensure a good representation of them, and not by criteria such as 
polynomial fits to the measured angular distributions, which are often truncated where the 
highest partial wave is still quite large. This is an important difference between our method 
and others 0J^,[rj,[n| which use observables in the form of Legendre series truncated at a 
point where the size of the highest coefficient is actually large. Our interpolation of the 
t-channel data used a method that makes no assumption about the number of partial waves 
which are important (see Appendix C of ref. ||), although we would expect higher partial 
waves to be progressively less important at lower momenta. 

At each £- value, sufficient terms were used to obtain a good fit to each amplitude and 
to ensure reasonable smoothness of the resulting helicity amplitudes from one t-value to 
the next. In practice, 4-5 J-values were used at the lowest momentum, rising to 6-7 at the 
highest. Initial estimates of Nj and Fj were obtained from || and used as starting values 
in a simultaneous fit to all the data at a given momentum. In Table | we show values of x 2 
per data point for the differential cross-sections and asymetries in the charged channel and 
for the amplitudes at each of the twenty momenta. At the higher momenta we also show 
values of x 2 P er data point for the neutral data. In Figs [I] and [| we show the fits and/or 
predictions for the 7r + 7r~ and 7r°7r° data at three representative momenta. To give some 
idea of the quality of the fits to the invariant amplitudes, we show in Fig. O the fits to the 
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real and imaginary parts of B and C as functions of t at a fixed value of the hyperbolic 
parameter, f = -0.04572 GeV 2 . 



IV. RESONANCES 

The extraction of "hard" information on resonances from partial-wave amplitudes is a 
notoriously difficult problem, which has not been solved in a rigorous way even for elastic 
two-body reactions, and ultimately one has to resort to plausible models. We will use criteria 
which have proved successful in analysing reactions such as ttN — > ttN. Thus, we will use 
combinations of Breit-Wigner resonances and flexible non-resonant backgrounds with initial 
parameters suggested by loops in Argand diagrams, enhancements in the integrated partial 
cross-sections, maxima in amplitude speed plots etc., not by fitting the original experimental 
data. Before doing that, however, it is useful to see to what extent the "continuity" of the 
partial-wave amplitudes, which were obtained by a series of fits at fixed values of the energy, 
is compatible with partial-wave analyticity. To test this we have fitted the "raw" amplitudes 
with the following parametric forms constructed by analogy with the work of ref . || . 

'iV_ N+ Nu 

n=0 n=l n=l 

where the variables z_ , z + and z u are given by 



Hj, L=J±1 = (l-z + ) L (l + z_ 



(11) 



z _ fl --(*-**) 1/2 (12) 



z _ a + -y4M 2 -t)V3 



and 



z = ^-^m>-ty/> 

u a u + (4m 2 -ty/ 2 - { V 

Here we use the values a_ = a u = 1 GeV and a + = 2 GeV for the mapping constants. In 



eq.[T| ti is given by 
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t L = m>(4-^- 2 ). (15) 

The three terms in eq.|ll] ensure that -fO,z,=j±i has the correct analyticity properties as a 
function of t. The factor (1 — z + ) L ensures the correct behaviour at the NN threshold while 
the factor (1 + zJ) 2 ensures a suitable high energy behaviour. We do not attempt to impose 
the correct behaviour at the itit pseudo-threshold since this is very far away from the energy 
region in which we are interested. The advantage in using such an a representation in the JL 
basis is that the NN threshold behaviour can be easily imposed, whereas in the helicity basis 
this is more difficult since the threshold behaviours of the two helicity states are correlated. 
The price for this decoupling is a correlation between the two JL basis amplitudes at t — 
to avoid a spurious singularity at this point. However, this is again so far from the energy 
region in which we are interested that we ignore this problem. 



The coefficients a n , b n and c n in eq.[Tl] are determined by fitting the values of i?j,L=j±i 
at the different t- values using the Pietarinen technique as described in ref. |§. In these fits 
we use AL = N u = 10 and N + = 15 and are able to achieve good fits; it should be remarked 
that in these fits we take no account of correlations, neither between different amplitudes 
nor between the real and imaginary parts in the same amplitude. The results of these fits are 
shown as the solid curves in Fig. f| for J < 5 where we show the dimensionless amplitudes 

hj±(W) = (pq) 1/2 H JiL=J±1 (W). (16) 

The J = 6 amplitudes are small and featureless and so are not shown. We conclude that, in 
the main, the single-energy amplitudes are compatible with partial-wave analyticity within 
plausible errors. In what follows, we have used these smooth amplitudes as the starting 
point for our extraction of resonance information. This is done for convenience, and it is 
worth remarking that no extra structure is introduced by the smoothing procedure and that 
we have checked that our conclusions are unaltered when the "raw" amplitudes are fitted. 

In Fig. ^ we also show the contributions of each partial-wave helicity amplitude to 
the integrated partial cross-section. Some systematic features are immediately apparent. 



Firstly, in a given energy region odd- J contributions are larger than those with even- J. This 
is of course directly related to the experimental observation that the 7r°7r° cross-section is 
approximately 1/3 of the tc~tc + cross-section over the whole of the energy range we consider. 
We also note the increasing importance of partial waves with lower L as the energy decreases. 
Many partial waves show counterclockwise loops, with the amplitude moving rapidly over 
some part of the energy range; the classic signal of resonance activity. To obtain resonance 
parameters, we have fitted the dimensionless amplitudes with the parametric forms 



where W = yt and 



2W-W min -W max 
x = . 18 

^max - W min 



In the background term the coefficients are complex parameters, and to ensure the 
correct behaviour at the NN threshold we set p = p/pb where pb is a suitable momentum, 
in practice taken to correspond to W = 2.1GeV. In the resonance term, the parameters are 
the mass Mr, the width T and the complex residue otj±. To ensure the correct behaviour 
at the NN threshold, we set 



«j± = 7J±(-) L+1/2 , P<Pr (19) 
Pr 

= 7J±, P > PR 

where pr is the value of p at W = M R and is a complex constant. Both amplitudes 
for a given J were fitted simultaneously, and in the absence of a realistic error analysis 
equal weights were assigned to each point. The results obtained by fitting the amplitudes as 
functions of energy, and the associated Argand diagrams, are shown as dashed lines in Fig. |] 
where it will be seen that the fits are as good as those using the more general parametrization 



of eq.ll 



The J = amplitude shows a clear resonance loop corresponding to a mass of 1.95 
GeV/c 2 and a width 0.16 GeV. The resonance dominates the real part at low energies, 

8 



but there is a large background at the upper end of the range. The imaginary part has 
significant contributions from both resonance and background terms throughout the entire 
energy range. In the case of the J = 1 amplitudes, there is a clear resonance signal in 
the L = amplitude corresponding to a mass of 1.97 GeV/c 2 and a width of 0.14 GeV. 
This is accompanied by a substantial background contribution. The L = 2 amplitude is 
by comparison much smaller. For J = 2, both h,2~ and /12+ amplitudes have resonance 
behaviour. The mass is 1.93 GeV/c 2 and the width is 0.15 GeV. The J = 3 state has a 
resonance in the h%- amplitude with a mass of 2.02 GeV/c 2 and a width of 0.22 GeV. The 
width in the case of the raw amplitudes is slightly larger at 0.26 GeV. The amplitude 
is very small. The resonance in the J = 4 wave, which again is more prominent in the 
L = J — 1 amplitude, has a mass of 2.02 GeV/c 2 and a width of 0.14 GeV, which rises 
to 0.26 GeV in the case of the raw amplitudes. The mass is in good agreement with that 
of the well-established f4(2050) resonance, and the values for the width span the accepted 
value of 0.21 GeV. The fact that this state emerges clearly from our analysis lends weight 
to the validity of our procedures and the parameters of the other predicted resonances. In 
the J = 5 wave there is evidence for a state at 2.19 GeV/c 2 with a width of 0.22 GeV, 
again coupling stronger to the L = J — 1 amplitude. Finally, there is some evidence for a 
broad resonance in J = 6, but as this is just outside the range of the analysis, and the raw 
amplitudes shown much more scatter, we do not comment on this possibility further and it 
is therefore not shown. 

The resonance masses and widths are summarised in Table Where two values are 
given, the first corresponds to fiting the raw amplitudes and the second to fiting the smoothed 
amplitudes. In the other cases both values are the same. Table || also gives the values of 
the product of the branching ratios 

Bj = B(R — > tt-r)B(R — > ~NN), (20) 

calculated from the values of the residue parameters jj±. The only case where there is a well- 
established resonance is the J p = 4 + state f4(2050) with a measured mass of 2.044±0.011 
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GeV/c 2 and a width of 0.208±0.013 GeV. If we identify our J = 4 state with this, we can 
use the known nn branching ratio of the f 4 (2050) of 17% to estimate the NN branching 
ratio to be between 2% and 8%, which is not inconsistent with a conventional qq resonance. 
If we assume that tttt branching ratios in the range 10-20% are also typical of daughter 
states in this mass range, then we can see from Table [Tl] that the other predicted resonances 
with J = 2, 3 and 5 have NN branching ratios similar to that of the f4(2050). However, the 
states with J = and 1 are different and have far larger branching ratios. Indeed for the 
J = state the coupling is probably unreasonably large, but in this case there is a very large 
background accompanying the resonance and in this situation our simple parameterisation 
may not be appropriate to extract accurate resonance parameters. Such states are unlikely 
to be conventional qq mesons, but are more likely to have a multiquark qqqq structure |L5] . 



V. OTHER ANALYSES AND CONCLUSIONS 



Three other analyses 1)12-1^] of the accurate tt tt + data [[□]] have been published. In ref. 



l~2|j it is assumed from the outset that the data can be fitted by partial-wave amplitudes 
each of which is parameterised in terms of towers of nearly degenerate resonances. The 
analysis then fits all the data simultaneously. The motivation for this assumption is 

based on qualitative features of the data, such as the persistence of a very large polarization 
over a wide energy range, but nevertheless is really only testing the compatibility of the 
data with an assumed parametric form, like analyses of earlier data QXJ] . The resulting speed 
plots show little evidence for resonance-like behaviour and only the J = 4 and 5 amplitudes 
show clear peaks. This illustrates the strong correlations that exist between the states in 
the various towers and makes it difficult to ascribe much reliability to the parameters of 
individual resonances. This is further illustrated by the predicted values of Bj which, with 
the sole exception of the J = 4 state, where the mass and width were fixed from experiment, 
all indicate very large branching ratios to the NN channel. 

The analysis of ref. |T3| is based on the Barrelet method M and therefore has the usual 
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problems intrinsic to that method associated with the unresolved phase. The initially very 
large number of possible solutions is reduced by using the 7r°7r° data and further reduced 
by imposing threshold behaviour on the way zeros appear, leaving just two solutions (A and 
B), both of which exhibit resonances which for (B) are in all waves from J = to J = 4. 

Finally, a simple amplitude analysis has been made [III] which finds no compelling evi- 
dence for any resonances. However, the analysis only used data in the restricted momentum 
range of 360 MeV/c to 988 MeV/c and is therefore not strictly comparable with the other 
analyses reported here. 



The work of ref. |T3[ is the most complete of the published analyses, but it is difficult 
to make a meaningful comparison with our own because of the significant difference in the 
input data to these two analyses i.e. our inclusion of invariant amplitudes and, via these, 
the requirement that the solution should be consistent with information in the itN — > ttN 
channel. To see the importance of this, we show in Fig. |3| the prediction of the solutions of 
ref. [[HJ for the invariant amplitudes as functions of t at the fixed value of £ = —0.04572 GeV 2 
together with our own fit to these amplitudes. It is clear, that despite the considerable errors 
on the invariant amplitudes, they do provide important constraints which are not satisfied 



by either of the solutions of ref. 



Nevertheless, there are some common results found in both analyses. Both find a clear 
signal for the f4(2050), and both find its coupling to the L = 3 amplitude is much stronger 
than to L = 5. Also, both analyses yield a low- lying J = state with a suggestion of 
abnormally large coupling to the NN channel. Agreement for other waves is much less good. 



For example, for J = 1 ref. ]T3[ finds a resonance only in solution B and it couples mainly 
to hi + , whereas the J = 1 resonance in our solution, although it has similar parameters, 
couples mainly to hi_. Likewise, there are substantial differences in masses, widths and 
couplings for the J = 2 and 3 waves. 

One place where it may be possible to decide between the various solutions is in the 
pp — > 7r°7r° channel where dcs data exist for only a limited momentum range. We show in 
Fig. ^| our predictions together with those of ref. |13|] for the differential cross section and 
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asymmetry in this channel at a lower momentum where there are no measurements. In the 
same figure we also show similar predictions at two higher momenta where measurements 
of the differential cross section exist [|3] . 

In conclusion, we summarise our overall findings. We have analysed data on pp — > titt, 
including very accurate data in the 7r~7r + channel, using constraints from crossing symmetry 
and analyticity to impose consistency of the solution with information on data in the irN — > 
ttN channel. Overall, we have found evidence for a rich spectrum of resonances whose 
couplings, with one or two exceptions, are in the range expected for normal daughter meson 
states with a qq structure. The exceptions are in low J states and have couplings far larger 
than expected. This may indicate the presence of more complex multi-quark structures with 



large couplings to the NN channel as expected in some models [15 
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FIGURES 

FIG. 1. Fits to the measured differential cross-section and asymmetry data for pp — ► 7r - 7r + . 



Our results are shown by solid lines and those of solution A of ref. [13] by dashed lines (solution B 
is very similar). (a) 404 MeV/c, (b) 1190 MeV/c, (c) 1500 MeV/c. 

FIG. 2. Predictions and/or fits to the measured differential cross-section and asymmetry data 
for pp — > 7r°7r°. Our results are shown by solid lines and those of ref. |]l3| by dashed lines, (a) 404 
MeV/c, (b) 1190 MeV/c, (c) 1500 MeV/c. 

FIG. 3. Fits to the real and imaginary parts of and as functions of t at a fixed value of 
the hyperbolic parameter £ = —0.04572 GeV 2 . The range of values for a given value of t from ref. 



|| are shown by the vertical lines, our fits are shown by solid lines and the predictions of ref. [13] 
by dashed lines. All values are in GeV units, (a) C(+), (b) bW, (c) CH, (d) BH. 

FIG. 4. Results of fitting the JL basis partial wave helicity amplitudes as functions of energy 
with parametric forms incorporating partial-wave analyticity based on eq.[ll] shown by the solid 
lines, and with simple resonance plus background forms, eq.|l^, shown by the dashed lines. The 
solid circles show the results of the single energy analyses, (a) J = 0, L = 1, (b) J = 1,L = 0, 
(c) J = 1,L = 2, (d) J = 2,L = 1, (e) J = 2,L = 3, (f) J = 3,L = 2, (g) J = 3,L = 4, (h) 
J = 4, L = 3, (i) J = 4, L = 5, (j) J = 5, L = 4, (k) J = 5, L = 6. 
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TABLES 

TABLE I. Values of \ 2 per data point obtained in fits to differential cross-section and asym- 
metry data for pp — ► ir~ir + , differential cross-section data for pp — ► 7r°7r° and amplitude data 



Momentum (MeV/c) -k tt + 7r°7r° Amplitudes 

360 1.89 0.03 

404 1.29 0.03 

467 1.76 0.05 

497 2.32 0.04 

523 1.68 0.05 

585 1.96 0.05 

679 1.70 0.04 

783 1.49 0.03 

886 1.05 0.03 

988 1.40 0.07 

1089 1.40 0.10 

1190 1.64 0.28 0.09 

1291 1.76 0.66 0.15 

1351 2.16 0.12 0.22 

1400 1.91 0.18 0.18 

1416 1.69 0.36 0.17 

1449 1.80 0.23 0.26 

1467 1.75 0.20 0.31 

1500 1.46 0.41 0.38 

1550 1.66 0.35 0.35 
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TABLE II. Resonance masses and widths in units of GeV, obtained from fitting partial- wave 
amplitudes as functions of energy, together with the values of the product of branching ratios 
Bj = B(R -► Tnt)B(R -► NN). 



J Mass Width Bj 

1.94, 1.95 0.16 0.18, 0.19 

1 1.96, 1.97 0.13, 0.14 0.053, 0.056 

2 1.93 0.15 0.013 

3 2.02 0.26, 0.22 0.013, 0.028 

4 2.02, 2.00 0.26, 0.14 0.013, 0.004 

5 2.19 0.22 0.004, 0.001 
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Figure 1(b) 
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Figure 1(c) 



2.00 — , 



CD 

E 
E 



0.00 — 




-2.00 



20.00 — 



a 

"D 
"D 



p, = 1 500 MeV/c 

Hab 



Experimental data from Ref. [11] 
Our fit 

Solution A of Ref. [13] 



0.00 




•1.00 



0.00 

cos e 



cm 



1.00 



Figure 2(a) 



1.00 —i 




Figure 2(b) 



1.00 — , 



CD 

E 
E 



in 



a 

"D 
"D 



0.00 



-1.00 



40.00 




20.00 — 



0.00 



p, =1190 MeV/c 

lab 



Interpolated data from Ref.[3] 
Our fit 

Predictions of solution A of Ref.[13] 




0.00 0.20 0.40 0.60 0.80 1.00 



cos e 



cm 



Figure 2(c) 



1.00 — , 



CD 

E 
E 



0.00 



-1.00 



10.00 



(f) 



a 

"D 
"D 



5.00 — 1 1 



0.00 



p, = 1 500 MeV/c 

Hab 



Interpolated data from Ref.[3] 
Our fit 

Predictions of solution A of Ref.[13] 




0.00 0.20 0.40 0.60 0.80 1 



.00 



cos e 



cm 



Figure 3(a) 




Figure 3(b) 




4.00 



5.00 



Figure 3(d) 



E 



10.00 



0.00 



-10.00 




10.00 



Our fit 

- + - Solution A of Ref. [13] 

— A— Solution B of Ref. [1 3] 



00 

CD 
DC 



0.00 



-10.00 — 







■ / 














7 


- 4 






4.00 



5.00 



Figure 4(a) 




-0.20 0.00 0.20 2.00 2.20 

Re E in GeV 

cm 



Figure 4(b) 




Figure 4(c) 





Figure 4(d) 




Figure 4(e) 




Figure 4(f) 





-0.10 0.00 0.10 2.00 2.20 

Re E in GeV 

cm 



Figure 4(g) 





Figure 4(h) 





Figure 4(i) 



0.05 — , 



0.05 — , 



0.00 



E 0.00 



-0.05 



-0.05 



-0.05 



0.00 

Re 



0.05 



2.00 2.20 

E in GeV 



cm 



2.00 



> 

CD 

CD 



LU 



E 

o 



2.20 




-0.05 



0.00 

Re 



100.00 — , 



£= 

-a 

"tr 

Q_ 



0.05 



50.00 — 



0.00 



J = 4, L = 5 




2.00 2.20 

E in GeV 



cm 



Figure 4(j) 





Figure 4(k) 



0.05 — , 



0.05 — , 



E 0.00 



O 



E 0.00 




-0.05 



-0.05 



-0.05 



0.00 

Re 



0.05 



2.00 2.20 

E in GeV 



cm 




